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Abstract 

We give a brief overview of recent results obtained through the gauge/gravity corre- 
spondence, concerning the propagation of a heavy quark in strongly-coupled confor- 
mal field theories (such as A/" = 4 super- Yang-Mills) , both at zero and finite temper- 
ature. In the vacuum, we discuss energy loss, radiation damping, signal propagation 
and radiation-induced fluctuations. In the presence of a thermal plasma, our empha- 
sis is on early-time energy loss, screening and quark-antiquark evolution after pair 
creation. Throughout, quark dynamics is seen to be efficiently encapsulated in the 
usual string worldsheet dynamics. 
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1 Introduction 

A substantial amount of evidence from heavy ion collision experiments at RHIC and, 
recently, LHC, supports the detection of the long-sought quark-gluon plasma (QGP), 
a hot and dense phase of deconfined strongly-interacting matter [1]. Energetic par- 
tons serve as important probes of this thermal medium, and there exists an enormous 
'jet-quenching' literature dedicated to analyzing the manner in which the plasma 
damps their motion, and is in turn disturbed by their passage [2]. Given the experi- 
mental indications that, in the range of accessible temperatures (of order a few times 
the deconfinement temperature), the QGP is strongly-coupled, perturbative QCD is 
believed to be inadequate for at least some aspects of the relevant calculations, cre- 
ating a demand for new theoretical tools. In the past five years, interesting steps 
have been taken towards meeting this demand via the AdS/CFT, or, more generally, 
gauge/gravity or holographic correspondence [31 H], starting with the seminal works 
[3 El m E] and continuing with a large body of work that has been reviewed in [9]. 

The gauge/gravity correspondence is by now well-established as a tool that grants 
access (often analytically) to a large class of strongly-coupled non-Abelian gauge 
theories, via a drastic and surprising rewriting in terms of string-theoretic (frequently 
just supergravity) degrees of freedom living on a curved higher-dimensional geometry. 
The gauge/gravity catalog includes some theories that display characteristic QCD- 
like physics such as confinement and chiral-symmetry breaking, but not, to date. 
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QCD itself. Thus, at our present stage of knowledge, contact with the real-world 
QGP produced at RHIC or LHC will be feasible only if QCD can be reasonably well 
approximated by at least one of the gauge theories whose dual description is known. 
In the past few years, encouraging signs in this direction have emerged even for the 
most rudimentary example [3], SU{Nc) maximally-supersymmetric (A/" = 4) Yang- 
Mills theory (MSYM), a conformal field theory (CFT) that at zero temperature is 
completely unlike QCD (and in particular, has a coupling that does not run), but at 
finite temperature is in various respects analogous to deconfined QCD. Besides the 
gauge field, MSYM contains 6 real scalar fields and 4 Weyl fermions, all in the adjoint 
representation of the gauge group. 

Following the very exciting findings on the shear viscosity in holographic thermal 
plasmas (which in turn led to calculation of many other transport coefficients 
prp. analyses of energy loss have been conducted considering various types of partonic 
probes of the MSYM plasma, including quarks O El [71 [12], mesons [131 [El [HI ITS] , 
baryons [T71 UHl [19], gluons [T71 [20] /c-quarks [T7| and various types of defects [2T] . 
Energy loss studies have also been carried out via the gauge/gravity correspondence 
in many other theories, including some with properties more akin to those of QCD — 
see, e.g., [221 [23] and references therein^. Other examples of drag force calculations 
may be found in the reviews [9]. While we are probably still far from achieving firm 
contact between experimental data and first-principles gauge/gravity calculations, 
the correspondence has already succeeded in enhancing our general intuition on the 
behavior of gauge theories at strong coupling, and has additionally provided useful 
suggestions for phenomenological models of the QGP (see, e.g., [251 [SEl [271 [23 [29]). 
In the present paper, we will review a number of recent results related to the evolution 
of a quark in a strongly-coupled gluonic plasma. 



2 Quarks, with Strings Attached 
2.1 Basic setup 

For our analysis we need two ingredients: a strongly-coupled thermal non-Abelian 
plasma, and a quark that traverses it. Via the gauge/gravity correspondence, a plasma 
in a (i-dimensional CFT is known to be described in dual language by a Schwarzschild 
black hole (more precisely, black brane) in {d+ l)-dimensional asymptotically anti-de 
Sitter (AdS) space. For concreteness, we will frame our discussion in the context of 
SU{Nc) MSYM in 3 + 1 dimensions, with coupling qym and temperature T, which is 
equivalent to Type IIB string theory on the (AdS-Schwarzschild)5 x geometry 



ds^ = G.mndx'^dx'' = — -hdf + rff^ + — + R^dQ^ , (1) 



^ -R^ o , 1 



1 4 ' IT ~ 9ym^c = a , Zh- 

zt It vrT 



^See also [24] for an interesting comparison between lattice results and its holographic counterpart. 
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(with a constant dilaton and units of Ramond-Ramond five-form flux through the 
five-sphere), where Ig denotes the string length, and Zh is the location of the event 
horizon. The radial direction z is mapped holographically into a variable length 
scale in the gauge theory, in such a way that z and z —> oo are respectively 
the ultraviolet and infrared limits. The directions = {t, x) are parallel to the AdS 
boundary z = and are directly identified with the gauge theory directions. The five- 
sphere coordinates are associated with the global SU{4) internal (R-) symmetry of 
MSYM. They will play no role in our discussion, so all results below will hold equally 
well in the more general case where the is replaced by a different compact five- 
dimensional space X5, which corresponds to replacing Af = A SYM with a different 
(3 + l)-dimensional CFT. 

Notice that the field theory temperature is identified with the Hawking tempera- 
ture Th = 1/ 'n'Zh of the black hole. In the limit of vanishing temperature {zh — >■ 00), we 
are left in ([1]) with a pure AdS geometry, which is dual to the (symmetry-preserving) 
vacuum of MSYM. The closed string sector describing (small or large) fluctuations on 
top of AdS fully captures the nonperturbative gluonic (-1- adjoint scalar and fermionic) 
physics. The generic (AdS-Schwarzschild)5 x geometry ([T]) is a special case of such 
a fluctuation, and is dual to a thermal ensemble in MSYM. The string theory descrip- 
tion is under calculational control only for small string coupling and low curvatures, 
which translates into iVc ^ 1, A ^ 1, i.e., a large number of colors and strong 
('t Hooft) coupling. 

It is also known that one can add Nf flavors of matter in the fundamental rep- 
resentation of the SU {Nc) gauge group by introducing in the string theory setup an 
open string sector associated with a stack of Nf D7-branes [SU]. We will refer to 
these degrees of freedom as 'quarks,' even though, being M = 2 supersymmetric, 
they include both spin 1/2 and spin fields. For \Nf <^ Nc, we are allowed to ne- 
glect the backreaction of the D7-branes on the geometry [221 131]; in the gauge theory 
this corresponds to working in a 'quenched' approximation that ignores quark loops. 
These branes cover the four gauge theory directions x'^, and are spread along the 
radial AdS direction from z = to z = Zm- The D7-brane parameter Zm is related to 
the Lagrangian mass m ^ \/XT of the quark through [5] 



1 2-iTm 



which simplifies to 



^oi(^ 



12> 



8 I 2m i 128 I 2m j 



(2) 



... . (3) 
27rm 

at vanishing temperature. For applications of this formalism to phenomenology, we 
must choose values of the mass parameter Zm based on the charm and bottom quark 
masses, m ^ 1.4,4.8 GeV. Following |32], at the temperatures relevant to RHIC this 
translates into z^/zh ~ 0.16 — 0.40 for charm and z^/zh ~ 0.046 — 0.11 for bottom. 
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An isolated heavy quark corresponds to a string extending up from the horizon at 
z = Zh to a. location z = Zm where it ends on the stack of Nf D7-branes. The string 
dynamics is governed by the Nambu-Goto action 

1 



Sng = J c^^o- V^iet^ = j (faCNG , (4) 

where Qah = daX^^d^X'^GMNi^) (a, ^ = 0,1) denotes the induced metric on the 
worldsheet. We can exert an external force F on the string endpoint by turning on 
an electric field F^i = Fi on the D7-branes. This amounts to adding to (jlj) the usual 
minimal coupling, which in terms of the endpoint /quark worldline x^(r) = X^(r, Zm) 
reads 

Sf = J drA.ixir))^^^ . (5) 

Variation of Sng + Sf implies the standard Nambu-Goto equation of motion for all 
interior points of the string, plus the boundary condition 

n^(r)|,=,„, = J-^(r) Vr, (6) 

where 11^ = OCng/ d{dzX'^) is the worldsheet current associated with spacetime 
momentum, and J^^ = F^yd^x'^ = {—jF ■ v,'yF) the Lorentz four-force. 

Notice that the string is being described (as is customary) in first-quantized lan- 
guage, and, as long as it is sufficiently heavy, we are allowed to treat it semiclassically. 
In gauge theory language, then, we are coupling a first-quantized quark to the glu- 
onic (-|- other MSYM) field(s), and then carrying out the full path integral over the 
strongly-coupled field(s) (the result of which is codified by the AdS spacetime), but, 
for the time being, treating the path integral over the quark trajectory x^{t) in a 
saddle-point approximation. (The effect of quantum fluctuations about the classical 
string conflguration will be discussed in Sections 12.21 and [3^ ) 

A crucial property of the gauge/gravity dictionary is that it identifles the end- 
point (s) of the string as being dual to the quark (or antiquark), while the body of the 
string codifles the proflle of the (near and radiation) gluonic (-1- other MSYM) fleld(s) 
set up by the quark. In other words, this dictionary teaches us that the usual 'QCD' 
string exists even for non-conflning theories, is infinitesimally thin, and actually lives 
in a curved 5 (+ 5)-dimensional geometry. The gluonic profile can be mapped out 
exphcitly by computing one-point functions of local operators {{Ti F'^{x)),{T^y{x)), 
. . .) in the presence of the quark, which, via the GKPW recipe [33], requires a de- 
termination of the near-boundary profile of the closed string fields {(f>{x, z), h^u{x, z), 
. . .) generated by the macroscopic string. 

It is important to keep in mind that the quark described by this string is not 'bare' 
but 'composite' or 'dressed'. In an interacting field theory, a particle of finite mass can 
only be regarded as strictly pointlike at zeroth-order in a perturbative description. 
In the nonperturbative framework made available to us by the gauge/gravity duality, 
a quark with finite mass automatically acquires a finite size. This can be seen most 
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clearly by working out the vacuum expectation value of the gluonic field surrounding 
a static quark located at the origin [M]. For m — )■ oo {zm — ^ 0), this is just the 
Coulombic field expected (by conformal invariance) for a pointlike charge. For finite 
m, the profile is still Coulombic far away from the origin, but in fact becomes non- 
singular at the location of the quark. The characteristic thickness of the implied 
non-Abelian charge distribution is precisely the length scale Zm defined in ([2]). This 
is then the size of the gluonic cloud that surrounds the quark, or in other words, the 
analog of the Compton wavelength for our non-Abelian source. 

Below we will also be interested in the description of a quark- ant iquark pair. The 
IIB strings we have introduced above are oriented, and a state with two oppositely 
oriented purely radial strings would correspond to a quark and antiquark that are 
merely superposed. With such boundary conditions, however, the configuration with 
lowest energy is given by a single U-shaped string that has both of its endpoints at 
z = Zm- This can again be verified by computing the corresponding gluonic profile in 
vacuum [351 ES], to verify that the falloff is more rapid than Coulombic, as expected 
for an overall color-neutral source. 



2.2 Late-time energy loss, Brownian motion and limiting ve- 
locity 

Combining the two preceding ingredients, we know that a heavy quark in a strongly- 
coupled MSYM plasma is described by a string on the Schwarzschild-AdS geometry 
([1]). As the string endpoint moves, its body lags behind it, exerting a drag on the tip 
that is the gravity-side realization of the damping force exerted by the plasma on the 
quark, as studied in OE]. For a restricted type of quark trajectories, this mechanism 
leads to energy loss at a rate 
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dt 2 



^T'^=^ ■ (7) 



V 



Various works have explored the way in which these heavy quark results are modified 
for heavy sources of the gluonic field in color representations other than the funda- 
mental (including the adjoint) [17], as well as for light quarks and gluons [371 138]. 

In the context of weakly-coupled QCD, it is known that energy loss at high par- 
ton velocity is dominated by the strong-interaction analog of bremsstrahlung, i.e., 
medium- induced gluonic radiation, while for low velocities coUisional loss becomes im- 
portant [2]. In the strongly-coupled MSYM setup made available to us by AdS/CFT, 
the perturbative terminology is no longer adequate. The fiow of energy /momentum 
along the body of the trailing string corresponds to energy/momentum transported 
away from the quark via the gluonic field generated jointly by the quark and the 
plasma, in a pattern that has been meticulously studied in a large body of work 
that began with [39], has been reviewed in [ID] and includes the recent additions 
[m Il2l SSI SH HSl 116]. For ease of language, throughout this paper we will continue 
to speak of radiation, even though this concept is not really appropriate within an 
infinite thermal medium. 
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We should stress that, within the framework of O [6], adopted also in this pa- 
per, the entire energy loss calculation is carried out inside the MSYM theory, whose 
coupling does not run and is taken to be large. For quarks that are extremely e- 
nergetic, the asymptotic freedom of QCD would lead one to expect deviations be- 
tween this scenario and the real-world QGP, and seek a hybrid approach describ- 
ing hard/perturbative emission of a gluon that subsequently propagates through a 
strongly-coupled medium, as advocated in [H1I1Z|I§ It is debatable, however, which of 
these two scenarios is more appropriate at the not extremely relativistic heavy quark 
energies achieved at RHIC. At the very least, the calculations discussed in this paper 
have a direct interpretation in terms of energy loss in a thermal plasma of strongly- 
coupled gluons and exotic matter ('XGP') of MSYM (or other CFTs), which is an 
interesting theoretical question in its own right. 

In O E] , the rate ([7]) was deduced for an isolated quark in the stationary or late- 
time regimes, i.e., assuming that the quark either moves at constant velocity, under 
the influence of an external force which precisely cancels the drag force exerted by the 
plasma, or has decelerated solely under the action of this drag force for a long period 
of time, and is about to come to rest. In the actual experimental setup, however, the 
configuration is neither stationary nor asymptotic: the quark is not externally forced, 
and slows down under the influence of the plasma, whose finite spatial and temporal 
extent imply that the late-time regime is not generically accessible]! Moreover, the 
real-world quark is not isolated, but is created within the plasma together with an 
antiquark. In Section|4]we will explore to what extent the rate of energy loss is affected 
by these issues. Also, as we will seen in Section [3l the above trailing string mechanism 
is in fact much more general, and accounts even for the radiative damping expected in 
vacuum [5ll[55]. That is, irrespective of whether a spacetime black hole is present or 
not, the body of the string plays the role of an energy sink, as befits its identification 
as the embodiment of the gluonic degrees of freedom. On the other hand, energy 
loss via the string does turn out to be closely associated with the appearance of a 
worldsheet horizon, as noticed initially in [ISl IH] at finite temperature and in [JH] 
(see also [STJ |58]) for the zero temperature case. 

The appearance of a black hole on the string worldsheet is crucial in particular 
to reproduce the expected Brownian motion of the dual quark in the hot medium. 
This connection was first worked out in detail in [5^ 1^ . for a static string on the 
Schwarzschild-AdS(i+i geometry, which as we know from Section 12.11 is dual to a 
static quark in a thermal bath of the CFT. At finite A, one must take into account 
small perturbations about the average string embedding. These are described by 

^The approach of [3 |47] involves the extraction from a hghthke Wilson loop of a transport 
('jet-quenching') parameter q, which characterizes the medium and whose value at strong coupling 
curiously differs from the customary definition of q as the average transverse momentum picked up 
per unit distance traveled 08] |49j . Some doubts about the procedure proposed in ^ i47| to calculate 
q on the gravity side of the correspondence had been raised by, e.g., [iSl [50l [51], but these concerns 
appear to have been resolved by a recent correction to the procedure [52], which happens to leave 
the result unchanged. 

•^This latter point has been emphasized from the phcnomenological perspective in [53] . 
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free scalar fields propagating on the worldsheet black hole geometry, and upon quan- 
tization, they are excited by Hawking radiation. This leads to Brownian motion 
of the endpoint/quark, whose detailed form is captured by a generalized Langevin 
equation. The authors of [59] reached this conclusion in arbitrary dimension by as- 
suming that the state of the quantized embedding fields is the usual Hartle-Hawking 
(or Kruskal) vacuum, which describes the black hole in equilibrium with its own 
thermal radiation. The authors of pD] focused on the case d = 4 and followed a dif- 
ferent but equivalent route, employing the dual relation between the Kruskal exten- 
sion of the Schwarzschild-AdS geometry and the CFT Schwinger-Keldysh formalism 
[6T| [62| [7] , together with the known connection between the latter and the generalized 
Langevin equation. These calculations were later generalized and elaborated on in 

[631 EH ESI EE]. 

Another very interesting prediction of AdS/CFT in this thermal context is the 
existence of a subluminal limiting velocity 




for the quark traversing the plasma. This follows simply from the fact that the 
quark velocity is dual to the coordinate velocity of the string endpoint, and v = Vm 
corresponds to a proper velocity at z = Zm equal to that of light [67]. A more general 
bound involving the external force F can be derived by requiring the Nambu-Goto 
square root to remain real [5B]. This same restriction on the velocity can be seen 
to arise from microscopic calculations of the meson spectrum [681 El]. Its validity 
for isolated quarks has been emphasized in [511 [56]. The existence of this limiting 
velocity might have interesting phenomenological consequences, such as the photon 
peak predicted in [27] or the Cherenkov emission of mesons analyzed in [25] . 



3 Isolated Quark at Zero Temperature 
3.1 Energy loss in vacuum 

To orient ourselves, we will first inquire into the rate of energy loss for an accelera- 
ting quark in vacuum, whose dual description involves the string moving on pure AdS 
spacetime. In this case we will have full analytic control over the system, which will 
allow us to develop some intuition on the problem at hand. Besides being of theo- 
retical interest in its own right, our vacuum analysis will provide a useful benchmark 
against which the finite-temperature results can be compared. 

A quark that accelerates in vacuum would be expected to emit chromoelectro- 
magnetic radiation. The first definite characterization of the radiation rate off an 
accelerating quark by means of AdS / CFT was worked out in an important paper by 
Mikhailov [70] . Remarkably, this author was able to solve the highly nonlinear equa- 
tion of motion for a string on AdSs that follows from (jl]), for an arbitrary timelike 
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trajectory of the string endpoint dual to an infinitely massive quark. In terms of the 
coordinates used in ([T]) (where for now h = 1), his solution is 

X'^{t,z) = z^^ + x^{t) , (9) 

with n = 0,1, 2, 3, and x^{t) the worldline of the string endpoint at the AdS boundary— 
or, equivalently, the worldline of the dual quark — parametrized by the proper time 
r defined through rj^y{dx^ / dT){dx'' / dr) = —1. 

From the structure of (jH]) we see that the behavior of the string segment located at 
a given time t and radial depth z (which according to our discussion in the previous 
section codifies the behavior of the gluonic field at the length scale z [71]) is completely 
determined by the behavior of the quark/string endpoint at the earlier, retarded time 
tietityz) defined by 

^ = ^^r^7=f + ^ret , (10) 
VI - f(trct)^ 

where the quark worldline has been parametrized by x°(r) instead of r, and v = 
dx/dx^. The definition ffTOj) can be shown to imply that tret is obtained by projecting 
back to the AdS boundary along a (fixed r) curve that is null on the string worldsheet 
[70] . in analogy with the Lienard-Wiechert story in classical electrodynamics. The 
embedding ([9]) thus describes a wave on the string that is purely outgoing: it is 
generated at the string endpoint and then descends along the body of the string, 
moving into the AdS bulk. This is why, among the infinite number of extremal 
string embeddings that satisfy the boundary condition associated with the given quark 
trajectory, the profile (Q is the one that is of physical interest for us: it is dual to a 
configuration where waves in the gluonic field move out from the quark to infinity. 

Working in the static gauge = t, = z, Mikhailov was able to reexpress the 
total energy of the string embedding ([9]) (via a change of integration variable z — )■ tret) 
as a local functional of the quark trajectory, 

m = ^ / dtrct— + E,{v{t)) , (11) 

where of course a = dv/dx^. The second term in the above equation arises from a 
total derivative that was not explicitly written down by Mikhailov, but can easily be 
worked out to be 



Vx ( 1 



2m=0 



. , /^^,/ 7m, (12) 

which gives the expected Lorentz-invariant dispersion relation for the quark. The 
energy split achieved in ( ITTl) therefore admits a clear and pleasant physical interpre- 
tation: Eq (associated only with information of the string endpoint) is the intrinsic 
energy of the quark at time t, and the integral over tret (associated with the body of 
the string) encodes the accumulated energy lost by the quark to its gluonic field over 
all times prior to tS Completely analogous statements can be derived for the spatial 



*This energy split has been extended to A:-quarks in the recent work [72] 
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momentum. No less remarkable is the fact that the rate of energy loss for the quark in 
this strongly-coupled non-Abelian theory is found to be in precise agreement with the 
standard Lienard formula from classical electrodynamics! The AdS/CFT correspon- 
dence thus teaches us that, in this very unfamiliar nonlinear setting, the energy loss 
turns out to depend only locally on the quark worldline. This feature has been argued 
in [29] to lead to an upper bound for the energy of a quark at finite temperature. 
In the next subsection we will see that both the rate of radiation and the dispersion 
relation are modified when the quark has a finite mass, and is therefore not pointlike. 

3.2 Radiation damping 

The radiation emitted by an accelerated charge inevitably backreacts on the charge. 
One effect, present already at the classical level, is a reaction force on the charge, 
that tends to damp its motion. And if the system is quantized, one additionally 
expects the emission of radiation to induce stochastic fluctuations of the charge's 
trajectory. The discovery of the gauge/gravity duality has opened the possibility of 
exploring these effects in the previously uncharted terrain of strongly-coupled non- 
Abelian gauge theories. We will discuss here radiation damping, and postpone the 
analysis of fluctuations until Section 13.41 

In the context of classical electrodynamics, and for an electron modeled as a 
vanishingly small charge distribution, the effect of radiation damping is incorporated 
in the classic (Abraham-)Lorentz-Dirac equation [73] , 



characteristic timescale set by the classical electron radius. The second term within 
the square brackets is the negative of the rate at which 4-momentum is carried away 
from the charge by radiation (as given by the covariant Lienard formula), so it is 
only this term that can properly be called radiation reaction. The first term within 
the square brackets, usually called the Schott term, is known to arise from the effect 
of the charge's 'near' (as opposed to radiation) field [711 C5]- The appearance of a 
third-order term in (fT3l) leads to unphysical behavior, including pre-accelerating and 
self-accelerating (or 'runaway') solutions. These deficiencies are known to originate 
from the assumption that the charge is pointlike. 

Analysis of radiation damping in the quantum non-Abelian case using traditional 
methods would be a serious challenge, but the AdS/CFT correspondence allows us 
to address it rather easily in certain strongly-coupled non-Abelian gauge theories 
[M| |55] . By Lorentz invariance, given the description of the static quark, we know 
that a quark moving at constant velocity corresponds to a purely radial string, moving 
as a rigid vertical rod. When the quark accelerates, the body of the string will trail 
behind the endpoint, and will therefore exert a force on the latter. Remembering 
that the body of the string codifies the SYM fields sourced by the quark, we know 




(13) 



with r the proper time, J^^ 



7(F • -u/c, F) the external 4-force, and te = 2e^/3mc'^ a 
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that in the gauge theory this force is interpreted as the backreaction of the gluonic 
field on the quark. In other words, in the AdS/CFT context the quark always has a 
'tail', and it is this tail that is responsible for the damping effect we are after. This 
is of course the same mechanism that yields the drag force exerted on the quark by 
a thermal plasma [5l [6], reviewed in the previous section. The analysis of [SU [55] 
makes it clear that the damping effect is equally present in the gauge theory vacuum. 

To obtain a noticeable damping effect, we must extend the analysis of the pre- 
ceding subsection to the more interesting case with Zm > 0, where the quark has 
a finite mass. As we have emphasized in Section 12. in this case our non-Abelian 
source is no longer pointlike but has size Zm- On the string theory side, the string 
endpoint is now at z = Zm, and we must again require it to follow the given quark 
trajectory, x^(r). As before, this condition by itself does not pick out a unique string 
embedding. Just like we discussed for the infinitely massive case, we additionally 
require the solution to be retarded, in order to focus on the gluonic field causally set 
up by the quark. As in [56], we can inherit this structure by truncating a suitably 
selected retarded Mikhailov solution. The embeddings of interest to us can thus be 
regarded as the z > Zm portions of the solutions In [SU |S5] it was shown that, 
with this information, the standard boundary condition for the string endpoint 
can be rewritten in the form of an equation of motion for the dressed quark, 

dr 2TTm. \ 27rm^ dr (AA}i 

~ 1 A T-2 ■ ^ ^ 




Notice that the characteristic length scale appearing in ( |T^ is precisely Zm = 
y/X/27rm, which as previously discussed, plays the role of the quark Compton wave- 
length. Let us now examine the behavior of a quark that is sufficiently heavy, or 
is forced sufficiently softly, that the condition ^ X\J^'^\/27rm? <^ 1 holds. It is then 
natural to expand the equation of motion in a power series in this small parame- 
ter. To zeroth order in this expansion, we correctly reproduce the pointlike result 
md'^x^ = J^^. If we instead keep terms up to first order, we find 

dr I dr 27rm? I 2Trm? dr 

In the 0{y/\) terms it is consistent, to this order, to replace J-"^ with its zeroth order 
value, thereby obtaining 

^yd^~2^~d^J~ ^^ll^~d^~d^' ^ ' 

Interestingly, this coincides exactly with the Lorentz-Dirac equation f|T3|) . with the 
Compton wavelength ([2]) playing the role of characteristic size for the composite 
quark. This is indeed the natural quantum scale of the problem. The radiation 
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reaction force in ( !T5|) is correctly given by the covariant Lienard formula, as expected 
from the result fill I) |7Uj . which we see then arising as the pointlike limit of the 
full radiation rate encoded in the right-hand side of fll4p . The Schott term in (1151) 
(associated with the near field of the quark), originated from the terms inside the r- 
derivative in the left-hand side of f|T^ . which we understand then to codify a modified 
dispersion relation for our composite quark. 

To second order in ^Jy\^\|2■Kvn? , we similarly obtain 




dr'^ dr^ dr 




d'^x^ d^x^ d^x" 
dr"^ dr"^ dr"^ 



7^ 



For compactness, we have grouped together all terms arising at the same order in the 
expansion, and used underbraces to mark the radiation reaction terms originating 
from the right-hand side of (IT^ (which now include corrections beyond the standard 
Lienard formula), and overbraces to indicate the near-field terms arising from the 
left-hand side of f|T^ (which incorporate corrections to the standard Schott term). 

One can continue this expansion procedure to arbitrarily high order in the ratio 
V^A|J^2|/27rm2. The full equation ([T4]) is thus recognized as a compact (reduced- 
order) rewriting of an infinite-derivative extension of the Lorentz-Dirac equation that 
automatically incorporates the of our non-pointlike non-Abelian source. It is 

curious to note that (1141) . which clearly incorporates the effect of radiation damping 
on the quark, has been obtained from ([9]), which does not include such damping for 
the string itself. The latter arises from the backreaction of the closed string fields set 
up by our macroscopic string, but these are of order l/N^, and therefore subleading 
at large N^. (though not necessarily at the real- world value A'^c = 3 [75]). 

The full physical content of f|T4|) can be made transparent by rewriting it in the 
form 

J"", (17) 



dP^ _ dp^^ ^ dP^ 



dr dr dr 

recognizing P^' as the total string (= quark -|- radiation) four-momentum, 

as the intrinsic momentum of the quark including the contribution of the near-field 
sourced by it (or, in quantum-mechanical language, of the gluonic cloud surrounding 
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the quark), and 




dr 



27rm2 




(19) 



as the rate at which momentum is carried away from the quark by chromo-electro- 
magnetic radiation. 

Unhke its classical electrodynamic counterpart ( IT3l) . the dressed quark equation of 
motion iHM has no self-accelerating (runaway) solutions: in the (continuous) absence 
of an external force, it uniquely predicts that the 4-acceleration of the quark must 
vanish. Interestingly, the converse to this last statement is not true: constant 4- 
velocity does not uniquely imply a vanishing force. 

It is natural to expect the energy split achieved in ( ITT]) or ( fT7|) to be somehow 
reflected in the geometry of the string worldsheet. It is interesting then that, on 
the string worldsheet dual to a quark undergoing arbitrary accelerated motion, there 
appears a black hole [56], whose event horizon naturally divides the worldsheet into 
two causally disconnected regions. The appearance of a worldsheet black hole had 
also been noted previously at finite temperature [IBl HH] , so the overall lesson is that 
such causal structure is intrinsically tied to energy dissipation, be it within a thermal 
plasma or in vacuum. It seems plausible then to interpret the regions outside and 
inside the black hole as corresponding respectively to the quark and the gluonic field, 
and postulate that the energy flow across this horizon should be related to the energy 
radiated by the quark [571 [56] . This possibility has been studied more closely in 
[SHI [77]. Unfortunately, this interpretation is not always correct: for general quark 
trajectories, one can show that the energy of the portions of the string outside and 
inside the dynamical horizon does not actually match the intrinsic and radiated energy 
of the quark [78] . 

3.3 Absence of broadening in the radiation pattern 

In the previous two subsections we have seen how a study of the string embedding 
yields important quantitative information on the dual quark dynamics, including its 
rate of energy loss, based on the natural split achieved in [70] of the total energy of the 
string. The latter is conserved on the fixed AdS (or AdS-Schwarzschild) background, 
but would of course decrease steadily if we take into account the gravitational (and 
dilatonic, etc.) radiation given off by the string in the course of its evolution. As we 
have noted, from the gravity perspective this radiation is suppressed by a factor of 
so it is consistent to neglect it at large Nc. On the other hand, through the 
GKPW recipe for correlation functions [33], it is this radiation (or, more precisely, 
the full metric perturbation produced by the string), evaluated at the AdS boundary, 
that determines the expectation value of the gauge theory energy-momentum tensor. 
This tensor contains not just the gross information about the total energy loss rate, 
but also the fine details about the directionality of the flow and the relative weight of 
the various dissipation channels. The problem of determining the spacetime profile 
of the gauge field disturbance produced by an arbitrary charge trajectory was solved 
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long ago by Lienard and Wiechert for classical electrodynamics, but, under various 
guises, remains of great interest today in the context of quantum non-Abelian gauge 
theories. It is therefore significant that the gauge/gravity duality has given us a 
useful handle on this problem for a varied class of gauge theories in the previously 
inaccessible regime of strong coupling. 

The translation between string and gauge theory disturbances was first explored 
in [33], which studied the dilaton waves given off by small fluctuations on an otherwise 
static, radial string in AdSs, and inferred from them the profile of the dual gluonic 
field observable (TrF^(x)) in the presence of an oscillating quark in vacuum. The 
MSYM waves were found to display significant temporal broadening, just as one 
would expect given that points on the non-Abelian field arbitrarily far from the quark 
can themselves reradiate. This feature emerges naturally in the gravity side of the 
correspondence, because motion of the string endpoint generates waves that move up 
along the body of the string, and each point on the string then emits a dilaton wave 
that travels back down to the observation point x on the AdS boundary, where, via 
the AdS/CFT recipe for correlation functions [33], the value of (TrF^(x)) is deduced 
by assembling together all such contributions, each with a different time delay. While 
informative, the results of [35j did not allow a definite identification of waves with the 
characteristic falloff associated with radiation, i.e., contributions that transport 

energy to infinity. The unambiguous detection of the latter calls for examination 
of the MSYM energy- momentum tensor, (T^,y(x)), which in the gravity side of the 
correspondence requires a determination of the gravitational waves emitted by the 
string [79] . 

In more recent years, motivated by potential contact with the phenomenology 
of the quark-gluon plasma [1], analyses of both (TrF^(a;)) and {T^u{x)) have been 
carried out in the case of a heavy quark moving at constant velocity through a thermal 
plasma (the only finite-temperature case where an exact solution for the corresponding 
string embedding is available [SIE]), in a large body of work that includes [HH] and has 
been reviewed in [ID]- The results are again compatible with the expected nonlinear 
dynamics of the (in this case, finite-temperature) MSYM medium. 

Given these antecedents, it came as a surprise when, back at zero temperature, 
additional calculations going beyond the linearized string approximation found, first 
for special cases [HI |12] and then for an arbitrary quark trajectory [H] (using the 
string embedding (^), that the ensuing energy density (Too(x)) displays no temporal 
broadening, and is in fact as sharply localized in spacetime as the corresponding 
classical profile. 

The apparent tension between the results of [35] and [121 HI] for the gluonic fields 
in vacuum was resolved recently in [IS], which carried out the (TrF^(x)) calculation 
beyond the linearized string approximation, and for an arbitrary quark trajectory, 
to put it on a par with the computation in [Hj. The results show that the gluonic 
near field at any given observation point depends only on dynamical data evaluated 
at a single retarded time along the quark trajectory (with luminal and subluminal 
propagation in the cases of infinite and finite quark mass, respectively). This proves 
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that, despite appearances, there is in fact no conflict between [35] and |121I13]. While 
it is true that the gluonic field emerges as a superposition of contributions with all 
possible time delays (as shown explicitly by the calculations in pT| W2\ HU US], and 
also by a recent reformulation of these as a super position of gravitational shock waves 
emitted by each point along the string |13]), it turns out that the net result evidences 
only the smallest of these delays. So, contrary to |12],|13], where it was suggested that 
the supergravity approximation is somehow leaving out the expected non-Abelian 
rescattering, in |1S] the no-broadening result was interpreted as a surprising prediction 
of the AdS/CFT correspondence for the net pattern of propagation in the CFT at 
strong coupling and with a large number of colors, under the assumption of a purely 
outgoing condition for the gluonic field generated by the quark. 

We should perhaps stress that this prediction refers only to propagation in the 
CFT vacuum. When a beam of radiation enters a plasma, it is expected to remain 
unbroadened only at distances that are small compared to the thermal wavelength, 
after which it will diffuse and eventually thermalize, consistent with the results of 



3.4 Fluctuations and the Unruh effect 

When going beyond the classical description of the string, two new effects are found, 
both of which are suppressed by a factor of the string length divided by the AdS 
curvature radius, or, equivalently (via (JT])), by an inverse factor of the CFT coupling 
A. On the one hand, we pick up the usual quantum fluctuations arising from the 
determinant of the path integral over string embeddings. These are present even for 
a static string (see, e.g., [80]), and lead to spontaneous deviations from the average 
endpoint /quark trajectory of the type studied, e.g., in [81]. On the other hand, the 
worldsheet black hole emits Hawking radiation, which populates the various modes 
of oscillation of the string. 

This second effect, which is present only for an accelerated trajectory and is thus 
associated with the quantum fluctuations induced by the gluonic radiation emitted by 
the quark, has been studied in [82] , for the special case of uniform proper acceleration 
A. In this case, the general solution (Q takes a form that was found independently in 
[35| 153] . The geometry induced on the worldsheet contains a black hole that is static, 
and thus amenable to explicit calculation. The analysis of the effect of Hawking 
radiation is in complete parallel with the story of Brownian motion [591 EHj recalled 
in Section \2.2\ except that it is complicated by the explicit time dependence present 
in the worldsheet geometry when presented in inertial coordinates. 

The problem simplifies if instead of working in the coordinates appropriate for an 
inertial observer we transform to a Rindler coordinate system adapted to an observer 
sitting on the quark. This transformation gives rise to an acceleration horizon both in 
the boundary and bulk descriptions. As a result, a state that is pure from the inertial 
perspective will generally be mixed from the point of view of the Rindler observers, 
because the field degrees of freedom accessible to the latter will be entangled with 
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degrees of freedom in the region beyond their horizon, which they must trace over. In 
particular, the pure AdS geometry expressed in Rindler coordinates, which is dual to 
the Minkowski vacuum of the CFT (as evidenced by the vanishing of the expectation 
value of the stress-energy tensor), is interpreted as a thermal bath at the temperature 
Tu = Ajlix. Contact is thus made with the celebrated Unruh effect [H], i.e., the 
fact that, from the point of view of a uniformly accelerated observer, the Minkowski 
vacuum behaves as a thermal medium. In Rindler coordinates, the string embedding 
is static and bends towards the Rindler horizon. It is well-known that the Rindler 
horizon of the CFT can be removed via a Weyl transformation, leading to the open 
Einstein universe. The corresponding bulk transformation drastically alters the radial 
foliation of the AdS geometry. Even though, by construction, in this new conformal 
frame the acceleration horizon is no longer visible in the boundary description, one 
can show that it is still present in the bulk, but lies at the fixed radial position that a- 
ccording to the AdS/CFT dictionary corresponds to the Unruh temperature. In other 
words, after this second transformation, the thermal character of the CFT state arises 
not from entanglement with degrees of freedom that lie beyond a spacetime horizon, 
but from the direct identification of the specific energy scale as the temperature of 
the CFT, in exact parallel with the dual interpretation of the Schwarzschild-AdS 
geometry. This exercise thus sheds light on the AdS implementation of the Unruh 
effect. (Related analyse s can be found in [581 ESI [85].) 

After the bulk diffeomorphism that removes the CFT horizon, the string embed- 
ding is not only static but also completely vertical. Interestingly, both the base string 
embedding and the background geometry at the location of the string are found to 
coincide exactly with the d = 2 thermal setup analyzed in [SH], which allows one 
to obtain the radiation-induced fluctuations of the quark simply by translating the 
results of that work. This close relation between the quantum fluctuations of the 
uniformly accelerated quark on Minkowski spacetime and the thermal fluctuations 
of a static quark in a thermal medium is evidently a direct consequence of the Un- 
ruh effect, but the reader should be aware that, for c? > 2, the detailed properties 
of this thermal medium are found to differ from those of the familiar homogeneous 
and isotropic thermal ensemble dual to the Schwarzschild-AdS^+i geometry. Back 
in the original inertial frame, the quark fluctuations are gove rned by Langevin type 
equations |H2]- 

4 Isolated Quark at Finite Temperature 
4.1 Constant velocity 

Having understood the rate of energy loss for a heavy quark that moves in the MSYM 
vacuum, in this section we restore z^, < oo — and consequently h < 1 — in the 
metric ([T]), to study the same quantity in the case where the quark moves through 
a thermal plasma. A thorough generalization of Mikhailov's analytic results [70] to 
this finite temperature setup would require finding the exact solution to the Nambu- 
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Goto equation of motion for the string on the AdS-Schwarzschild background, for 
any given trajectory of the string endpoint at > 0. Sadly, this has not yet been 
accomphshed. Nevertheless, based on the results discussed in the previous section, 
we expect the total energy of the string at any given time to again decompose into 
a surface term that encodes the intrinsic energy of the quark and an integrated local 
term that reflects the energy lost by the quark to the thermal medium. 

There are two easy cases where one can show analytically that this expectation 
is borne out the static quark, where there is of course no energy loss, and the 
quark moving at constant velocity. An important difference with respect to the T = 
case analyzed in the previous section is that here the surface contribution arises not 
only from the lower {z = Zm) but also from the upper {z = Zh) endpoint of the string. 
Since Zh marks the position of an event horizon, for any finite coordinate time t the 
value of the surface contribution at Zh is not influenced by the behavior of the z < Zh 
portion of the string, but depends only on the string's configuration at t — )■ —oo. 
The same interpretation can be then carried over to the gauge theory: for arbitrary 
finite-temperature configurations, the terms in the quark dispersion relation arising 
from the upper string endpoint will encode a contribution to the energy of the state 
that depends solely on the initial configuration of the quark+plasma system. 

4.2 Early-time energy loss 

As we mentioned in Section \2.2\ expression ([7]) was derived under a number of sim- 
plifying assumptions. The one that matters most for our purposes is the restriction 
made by the authors of [21 E] to a configuration where the quark is either moving 
with constant (possibly relativistic) velocity as a result of being pulled by an external 
force that precisely balances the drag, or is unforced but moving nonrelativistically 
and about to come to rest. Either of these scenarios requires the interaction be- 
tween the quark and the medium to occur over a considerable period of time. The 
actual energy loss might thus be expected to differ from ([7]) in a situation where the 
quark moving through the plasma is accelerating, or in the initial period following its 
production within the thermal medium. This point has been emphasized from the 
phenomenological perspective in [53]. The estimates there are based on perturba- 
tive calculations, so it is interesting to inquire into this issue in the strongly-coupled 
systems available to us through the AdS/CFT correspondence. 

Now, the restriction in [S] [6] to the stationary or asymptotic cases was of course 
implemented to gain analytic control on the problem of energy loss. The authors 
of [77] were able to slightly extend the stationary result to the case of a slowly de- 
celerating quark, working in an expansion in powers of \f\T jm. Away from these 
regimes, it is difficult to follow the evolution of the quark in the thermal plasma, or 
equivalently, of the string on the AdS black hole geometry. 

The first gauge/gravity exploration of early-time energy loss was carried out in 
[56] , by numerically integrating the Nambu-Goto equation of motion for the string 
dual to a quark initially at rest, which is accelerated along one dimension by a time- 
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dependent external force that is turned off after a short interval, allowing the quark to 
move thereafter only under the influence of the plasma. The results show a qualitative 
difference between the initial stage (0 < t < treiease) where the quark is accelerated by 
means of the external force F{t), and the second stage (treiease < ^ < ^breakdown) where 
it moves only under the influence of the plasma. In the former stage, the rate of energy 
loss, for values of m in the neighborhood of the charm or bottom masses, is in fact 
nearly identical to the corresponding rate in vacuum, given by the modified Lienard 
formula f[T^ . In the latter stage, which would appear to be more relevant from the 
phenomenological perspective, the quark was found to dissipate energy at a rate much 
lower than the late-time result ([7]). Indeed, for the values Zm/zh = 0.2,0.3,0.4 that 
(as noted below ([3])) are appropriate for the charm and bottom quarks, the late-time 
friction coefficient /i = —{l/pg){dpq/dt) is respectively fii^tc/'^T = 0.25,0.41,0.59, 
but the numerical results for v(t) at treiease < ^ < ^breakdown are best approximated 
by /ieariy/7rT = 0.08,0.15,0.26. A more detailed discussion of the energy lost by the 
quark in both stages was given in [56], and a similar comparison of vacuum versus 
thermal energy loss rates was performed later in [86], for the special case of a quark 
undergoing uniform circular motion. 

Additional control on the problem of early-time quark damping in a strongly- 
coupled plasma was gained recently in [87], which studied thermal effects analytically, 
as a small perturbation on the zero-temperature evolution. This is reasonable on 
general physical grounds during the initial stage of motion through the plasma, and 
is moreover supported by the numerical results of [56]. Interestingly, this framework 
naturally makes contact with the previously known limiting velocity Vm (defined in 
([H])) for a quark immersed in the strongly-coupled plasma. The results of [H7] identify 

1Tl.Vm 1 /T^ / s 

Eg = - -v^T (20) 

as the thermally corrected intrinsic energy of the quark, and 



dt 2(f^ — 16m2(w^ — v^) 



(21) 



as its rate of energy loss. 

For parameter values appropriate to RHIC, this early-time rate can be up to 
six times smaller than the late-time result ([7]). This comparison is consistent with 
the (more limited) numerical results of [56] described above. The functional form 
of the results is also of interest: from the derivation in [87] it becomes clear that 
the T^-dependence seen in ( 12T|) (and anticipated in an estimate based on saturation 
physics [57]) follows naturally from the first thermal corrections in the metric, and 
via dimensional analysis, this fixes the accompanying dependence on the quark mass, 
which in turn explains the peculiar power of the 't Hooft coupling. Another prominent 
feature of this formula is its velocity-dependence: whereas the late-time friction coeffi- 
cient deduced from ([7]) is momentum- independent, its early-time counterpart depends 
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strongly on momentum. This feature should constitute an interesting experimental 
signature when used as input for phenomenological models, as in [23 [22] . 

The rate (121 p is limited to the initial stage where the quark is only starting to feel 
the effects of the plasma, and should be reliable for times that are not too close to the 
lifetime of the real-world QGP. This suggests that the experimental results should 
effectively arise from some sort of average or interpolation between the early and late 
energy loss rates (ET]) and ([7]). For times of order XjixT^ v^^ — v^, one expects the 
string embedding to be significantly modified by the appearance of and higher 
corrections neglected in the linearized analysis of [87]. As the string continues to 
evolve, a dynamical worldsheet horizon (located z = in the remote past) is 
expected to move up toward smaller values of z^ and then (if the quark is undisturbed) 
move back down towards the spacetime horizon [56] . When the wavefront on the string 
crosses the worldsheet horizon, one would have to deviate from the purely outgoing 
condition used in [87], to accommodate the fact that, as the medium is disturbed 
by the quark, it is expected to radiate gluonic waves back towards the source. The 
net result should be to yield a damping coefficient /i(t) that smoothly interpolates 
between the early-time result and the late-time form. 

5 Pair Creation within the Plasma 
5.1 Quark-antiquark evolution 

As was mentioned in Section [2^ in the experimental setup the heavy quark is created 
within the plasma together with its corresponding antiquark, and the presence of the 
latter would be expected to substantially modify the gluonic fields in the vicinity 
of the quark, consequently affecting its evolution. The evolution of a heavy quark 
and antiquark that are created within the plasma and then separate back to back 
was first studied in [5], and explored in greater depth in [56] . In dual language, this 
corresponds to a string with both of its endpoints on the D7-branes at ^ = 2^, such 
that the endpoints are initially at the same spatial location but have initial velocities 
in opposite directions. 

In setting up the problem, one realizes that on the gauge theory side there are 
actually two distinct quark-antiquark configurations: the product of a fundamental 
q and an antifundamental q can lead to a q-q pair either in the singlet or the ad- 
joint representation of the SU (iVc) gauge group. There is a counterpart to this in 
the gravity side, because there are precisely two distinct types of consistent initial 
conditions for the string, which lead it to expand into either a U-shape or a V-shape 
as its endpoints separate [56]. Curiously, the initial quark velocity is found to be 
freely adjustable in the adjoint, but not the singlet, configuration — in the latter case 
it is invariably fixed at the limiting value v = Vm- 

As the quark and antiquark separate, we expect them to eventually enter the 
late-time regime where their rate of energy loss is given by ([7]). The agreement be- 
tween the analytic and late-time numeric results is most cleanly seen if instead of 
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comparing graphs of x(t) or v(t) for the quark (where one would need to look at 
t — )■ oo), one examines the plots of v{x). These plots are shown in Fig. [H for mass 
parameter Zm/zh = 0.2 (in the neighborhood of the charm quark), corresponding to 
a limiting velocity Vm = 0.9992. It is evident from the figure that the late-time be- 
havior is well-described by ([7]) in all cases, but there is also an initial period where 
the behavior is different. This difference is clearly more significant for the singlet 
than the adjoint case, just as one would expect. For values of the mass in the neigh- 
borhood of the charm or bottom masses, the numerical results of indicate that 
the initial evolution is essentially identical to what it would be in the absence of the 
plasma. This picture seems rather close to the phenomenological discussion given in 
[53] (in the context of collisional energy loss): when the singlet quark- ant iquark pair 
is formed within the plasma, there is a delay before the interaction between the newly 
created sources and the plasma can set up the long range gluonic field profile that 
is responsible for the late-time dissipation. To examine in more detail the transition 
to the late-time behavior, for a variety of trajectories one can determine the point 
{xf,Vf) beyond which the numeric v{x) curve agrees with the late-time curve to a 
given accuracy /. (The work [SS] considered / = 5 or 10%.) These transition points 
are schematically indicated by the arrows in Fig. [T] 
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Figure 1: Quark evolution (velocity as a function of traveled distance, in units of 
l/vrT) for five different initial conditions. The dotted curves show the results of the 
numerical integration in [56], contrasted against fits in solid red that use the late-time 
analytic expression ([7]). The three dotted curves starting at Vq = Vm. describe singlet 
configurations with different initial velocity profiles and energies. Notice that the 
purple curve describes a situation where the quark and antiquark turn around and 
come to rest while approaching one another. The two remaining curves arise from 
adjoint configurations with different energies and initial quark velocities. The vertical 
arrows mark the transition distance where each trajectory enters the late-time regime. 
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5.2 Screening length 

By the time when the quark moves beyond Xf{vf) and therefore enters the late- 
time regime, it is certainly insensitive to the presence of the antiquark. It is natural 
then to ask how the transition distance Xf{vf) compares against (half of) the length 
-^max(f ) beyond which the quark and antiquark are screened from each other by the 
plasma. Based on what we just said, we know that for a given velocity we must have 
-f'max/2 < Xf, but the actual comparison informs us on whether the transition to the 
regime where the quark experiences a constant drag coefficient occurs right after the 
quark and antiquark are screened from each other by the plasma, or if there is an 
intermediate regime where the quark moves independently from the antiquark but 
nevertheless feels a drag force that differs from the stationary result of El E] , as we 
found when studying the early-time evolution of an isolated quark in Section 14.21 

The screening length for infinitely massive quarks in MSYM was computed in 
[m [15], by considering a quark-antiquark pair moving jointly through the plasma (a 
related length was plotted in |T3], and [15] gave in addition an extensive analysis of 
the energy of the q-q pair). This entails a study of a U-shaped string whose endpoints 
move in the same direction along the Schwarzschild-AdS geometry, a configuration 
that, interestingly, feels no drag in the —> oo limit, on account of its being color- 
neutral [131 El [IS] El Over the entire range < v < 1 the screening length may be 
approximated as [15] 

w^o^^a-^Y^s (22) 

while in the ultra-relativistic limit, it can be shown analytically that [11] 

1 '^-3/44 3/2 n 740 

The 1/4 in the exponent here has an intuitive interpretation in terms of the boosted 
energy density of the thermal medium [1?]. The full curve L^i^^{v) does not deviate 
far from this asymptotic form, so a decent approximation to it is obtained by replacing 
0.743 — )■ 0.865 in (^^, to reproduce the correct value at f = (at the expense of 
introducing a 16% error as f — 1) |14] : 

^.ax(.)-^(l--^)^/^ (24) 

A comparison between the two approximations (!22|) and (!24|) is shown in Fig. |2) 
overall, the exponent 1/3 is better than 1/4 in the sense that it implies a smaller 
squared deviation from the numerical results, even though 1/4 leads to a smaller 
percentage error in the range v > 0.991 (7 > 7.3). An attempt to better parametrize 
the deviation away from the ultra-relativistic behavior was made in [SB]. In any case, 
one should bear in mind that the region of principal interest at RHIC is not really 
f — 7- 1, but 7f ~ 1. 



^The 0{1/N^) drag force has been worked out by indirect means in [16]. 
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Figure 2: Screening length Lmax (in units of l/27rT) as a function of velocity (in 
black) compared against the approximations ( 122|) (in red) and ( 124|) (in blue). 



To compare against the transition distance Xf defined in the previous subsection 
for the back-to-back quark and antiquark, the authors of [56] extended this calculation 
to the case of finite mass. The main difference with respect to the m — !■ oo case 
is that now the 'ultra-relativistic' region would refer to the limit where the pair 
velocity approaches the limiting velocity Vm < 1. The screening length Lmax('y) is still 
relatively well approximated in the full range < f < by the natural modification 
of the 2;^ = fit (!22|) (and in all cases, the fit analogous to ( 12^ does a poorer job). 
The behavior in the v ^ Vm region can still be determined analytically, and turns out 
to be linear in v'^—v'^. Over the full velocity interval, the screening length Lmax(f ) and 
the transition distance Xf{v) and (half of) are found to be of comparable magnitude 
and scale with velocity in a similar manner |56]. A similar agreement was reported in 
|38j . and later related work may be found in |86]. Notice that the agreement holds in 
spite of the fact that the two relevant string configurations are quite different, with 
motion respectively along and perpendicular to the plane in which the string extends. 

The conclusion then is that the transition to the late-time regime ([7]) takes place 
immediately after the quark and antiquark lose contact with one another. That 
is to say, unlike what we found for the forced isolated quark in Section 14. 2[ here 
there is no appreciable intermediate stage where the quark and antiquark decelerate 
independently from one another at a rate that differs substantially from the late-time 
result of O [6] . This then extends and at the same time delimits the region where the 
analytic results of O [6l [7] can justifiably be used to model energy loss in heavy ion 
collisions. 

The way in which the findings reviewed in the present section interface with those 
of Section 112] was discussed in [87]. The upshot is that, in the more realistic situation 
where the presence of the accompanying antiquark is taken into consideration, the 
early-time damping results would be most relevant in the stage before the separating 
quark and antiquark lose contact with one another. This is natural, because screening 
emerges only when thermal effects are substantial. For a quark ploughing through a 
significant portion of the plasma, then, the late-time formula ([7]) actually stands a 
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chance of controlling the larger portion of the evolution. 
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